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Abstract 

We show that the energy levels predicted by a 1/iV-expansion method for an N-dimensional 
Hydrogen atom in a spherical potential are always lower than the exact energy levels but 
monotonically converge towards their exact eigenstates for higher ordered corrections. The 
technique allows a systematic approach for quantum many body problems in a confined 
potential and explains the remarkable agreement of such approximate theories when compared 
with the exact numerical spectrum. 
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A fundamental theoretical problem in suggests that an increase in the number of 

the realm of many body physics concerns degrees of freedom often simplifies the theo- 

the technical difficulty in making precise retical analysis jl| . A perturbative approach 

theoretical evaluations of physical observa- requires at least one dimensionless parame- 

tions, even more so in problems involving ter and if we couple this fact with the pre- 

quantum systems. In most cases, this is vious statement, it effectively implies that 

compounded by the practical difficulty in es- as we go on increasing the dimensionality 

tablishing a suitable approximation scheme of this parameter, the perturbation analysis 

that conjoins simplicity with effectiveness, becomes more and more simple Of- 

An interesting observation in connection ten it is found that a problem of inherently 

quantum mechanical origin can be mapped 
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on to a classical phase space in the N — > oo 



1 



limit thereby reducing a quantum problem Although the model is nothing new jlfl| . 

to a classical one In other words, one however our objective here is. We intend to 

then has a limit where quantum interfer- study the efficacy of this expansion method 

ence effects simply die out paving the way by calculating the energy eigenvalues and 

for a simple classical analysis. The excited showing that to each order of correction, the 

states for such a system can be obtained large-N expansion method always predicts 

as an expansion in 1/N around the mini- a slightly lower potential as compared to 

mum of the effective classical potential V e g. the exact eigenvalue obtained numerically. 

Such an approach is not at all uncommon in This is remarkable since this implies a cer- 

statistical physics 0, 0| in problems which tain monotonicity in these perturbation cor- 

allow for at least a minimum. In many of rections which tells us that the corrections 

those cases, the large-N limit has been fruit- are always positive, a fact that has often 

fully utilized in dealing with equilibrium as been tacitly assumed in related calculations 

well as non-equilibrium problems in classi- lj| 12 1- We argue that this is the underly- 

cal critical phenomena [l]. In quantum me- ing reason which makes this method more 

chanics too, 1/N expansion method has a dynamic compared to standard perturba- 

long precedence. Detailed accounts of re- tion technique which is limited strictly to a 

lated applications can be obtained from re- weak- coupling regime. In a following work, 

view articles like the one due to Chatterjee we build on this principle and analytically 

as also from more informal nar- solve for the three-body problem of interact- 

ratives like ]} 0| • The versatility and flexi- ing electrons using an exact Coulomb poten- 

bility of this technique has allowed it to be tial [18^ . 
used in a range of diverse topics, starting 

c n u ,i ]••!•! In the first paper of the paper, we do a 
from held theoretic studies m high energy 

rehash of the N-dimensional quantum me- 



physics [1^ 



151 ] to problems on earthquake 



, 10 l , , . , chanics for a single electron in a spherical 
dynamics |lo|J as well as on problems m col- 
loidal physics ^| confining potential and then defining the 

potential in the relative frame of reference, 

In this brief report, we shall use the we go on to solve the stationary state prob- 

1 /iV-expansion method to study the prob- lem using the 1/iV-expansion method. As 

lem of an N-dimensional Hydrogen atom already stated, we then proceed to calculate 

confined in a Harmonic oscillator potential, the energy corrections due to this method 



for both ground and excited states and show 
that all higher ordered corrections have a 
steady monotonicity that ensures a large-N 
eigenstate below its exact (meaning experi- 
mental) counterpart. 

Taking cues from standard literature 8, 

n 

HO], we begin with the Hamiltonian for the 
center of mass of an N-dimensional electron 
in a spherical potential 



X\ = r cos#i sin 2 sin 3 ... sin#/v-i 
22 = rsin#isin#2sin#3...sin#/v_i 
x 3 = r sin 6*2 sin 6*3 sin 6* 4 ... sin A r_ 1 



Xfc = r cos8k — 1 sin 9k sin &k+i--- srn 6n-i 



H 



-n 

2m f 



+ V N (f 



(1) 



Using standardized units h = m e = 1 
(h=Plank's constant and m e =mass of the 
electron), the Hamiltonian can be rewritten 

as 



H 



1 



V N 2 + V N (r) 



where terms have their usual meaning. 
The potential being radial Vjv(r) = Vjv(r) 
and this gives the eigenvalue equation 



xn-i = r cos6 l Ar_2 sin 9n-i 
xn = tcos9n-i 



(4) 



where r is the r is the radial distance and 
9k(k < N — 1) are the angles defining the 
hyper-spherical space, 9^-i being the az- 
imuthal angle. ijj(r) is the eigenfunction of 
this system. The above definition can now 
be used to obtain the radial equation of mo- 
tion \m] 



H^[f) = [-- V N 2 + V N (r)}ij(r) = E^{r) 

(3) 

For a system with spherical symmetry, 
the curvilinear coordinates can be written 
as follows (generalisation of the treatment 
available in [191 ]) 



r 1, d 2 N — 1 d . 1(1 + N -2) 

| — ( 1 )+ 

2 dr 2 r dr 2r 2 

+ V N (r)]R(r) = ER(r) (5) 

where /'s are the angular quantum num- 
bers and R(r) is the radial wave func- 
tion. Using the transformation u(r) = 
r ( N ~ 1 )/ 2 R(r), we can now absorb the first 



derivative in eqn. (0). The reconstructed One can now add quantum fluctuations and 
radial equation of motion is now given by study the behavior of the system close to 



the classical minimum r 



We go be- 



1 d 2 R (1 — -)(1 — -) Vk(t) yond this description in the sense that we 
'071+^ r — ^ — rr~\ u { r ) = Eu(r) 

2 dr ' 8r k consider a finite sized electron instead of a 

(6) 

fixed mass and consider fluctuations around 



In the above, we have used k = N + 21. 
At this point, the meaning of the large-N 
limit turns out to be pretty obvious. It 
means that k — > oo (since N is large) en- 
compasses the idea of a stationarity limit 
for a very heavy classical particle of effec- 
tive mass k 2 where the particle is localized 
at the point r = r , the point r in turn 
defining the minimum of the classical po- 
tential V e ff = + ^2 ■ The ground state 
energy of such a localized system is given by 



the classical stable minimum. To do this we 
revoke the original radial equation eq. © 
prior to the large-N limit being imposed on 
it. Using the 1/N expansion technique, we 
now embark on a stepwise evaluation of the 
energy eigenvalues due to the quantum fluc- 
tuations close to the classical minimum. We 
define the eigenvalue problem as follows 

[H + V(r)}ij(r) = E^(r) (7) 



£ _ , r \ The ground state eigenvalue equation 

w • , .n r c ,i H ib(r) = Enibir) has already been de- 

We now consider a specific form for the rv ' urv ; J 

, , • i f , • tt i \ j j , , fined through equation Q (En = while 

potential function V N [r) and proceed to cal- ° H v 2/ 

i , ,i i • i j . ,i yfr) is the part of the Hamiltonian that 

culate the higher order corrections m the v ' 1 

i xt v •, rpi ii i f contributes to the quantum fluctuations, 

large- IN limit, the model we choose for H 

• .„ ,-,ii • Taylor's expansion allows this perturbation 

the purpose is an oscillator with anharmomc j v r 

a , ,. rr, r ,1- i • i Hamiltonian to be represented as 

fluctuations, lhe reason tor this choice has v 

been accentuated by the observation that 

such a description, albeit simple, yet is able y _ V"(y ) | (r ro)V'(ro) I — — V"(tq) I . 

to reproduce a good estimate for the en- ^ 

ergy eigenstates |20] when compared with where the pr i mes denote derivatives with 

numerical |2l| as well as with experimen- reS pect to r. Before proceeding any ur- 

tal [23 result. For a simple harmonic os- ther; we make a variable transformation 

dilator V N (r) = \uj 2 r 2 which gives r = from r _ x where x = vj( r _ r ^ and 

V eS = |^| and eventually E M = |^|. transform eq. (0) likewise. In the translated 



coordinate system, the complete eigenvalue Higher-ordered parameters like e 3 have 
equation is given by non-zero values for anharmonic oscillations. 

The above description allows us to re-frame 

1 d 2 u 4 3 x x 2 an effective classical potential V e s in the 

" ol~2 + M(l -T + 15X1 - 2^= + 3- 

^ ux K k \J k k large-N limit but now including higher- 

3 

— 4 — j- + ...) + r 2 {V(r ) + r V'(ro)(l H — ordered fluctuations. It has the form 
+ ! ^V r "(ro)(l + ^ + 2^)...}] M (x 



* W = 4t^? + >) (13) 

r 2 u(x) (9) 



^ a , x „/ cttV 2 k k k 

k 



where V(R) represents some oscillator 
In the analysis of the above equation we potential having a minimum at R ^ a point 

consider all terms up to 0(r 2 ) and evalu- which c&n be obtained from the relation 

ate coefficients for increasing powers of x 

starting with x°. A little rearranging now 9 . . . 

-Q^Ves{n)\ R=Ro = (14) 

allows us to rewrite the eigenvalue equation 

Defining the potential as in eq. (JTTj) and 



in terms of the variable X clS follows 

[H + V{x)]ip(x) = \ip{x) (10) 

where 



then applying the optimization criterion as 
in eq. (fTljl , we arrive at the quadratic equa- 
tion 

3e 3 i? 2 - Vkuj 2 R + e x = (15) 

Id 2 1 

H = H — cu 2 x 2 + 6q which gives the solution = 

2 doc 2 -y^fc ^ ^ /fe 4 

- , , 1 , o, 1 , n 4 . — - — ^-r— — ■ To check the stability at 

V(x) = —j=(eix + e 3 x s ) + T (e 2 x 2 + e 4 x 4 ) & ^ J 



the point R = Rq, we evaluate the second 



Vk k 

~fc3/2 ^ X + ) (1-0 derivatives and find that the two roots of eq. 



, > ,e\ 1 i r 1 • flTBT) satisfy the relation 

where A = (tJTo an d for a harmonic v y 

oscillator potential the constants e& and 5k 

are given by 



w 2 + 6e 3 , v 7 ^ 2 ± y/kuo 4 - 12ei 



8 2 8A; 64 ^he above result implies that the minima 

ei = 1, e 2 = -3/2, e 3 = -r V"(r ) -(l^ are subject to the restriction koo 4 > 12eie 3 . 



An idea of the exactitude of this analysis can 
be had from an evaluation of the parameters 
using a simple harmonic oscillator potential. 
This gives e% = —1/2, thereby naturally val- 
idating the restriction. The conclusion re- 
mains unchanged even after adding higher 
ordered anharmonic terms to the potential. 
To leading order in expansions, we now have 
the large-N expanded energy eigenvalue for 
the ground state as follows 

e = ^° (+)2 +4^ (+ W 3 ;?o (+)2 )+? 

2 ro ?"o k 

(17) 

where r 2 = j-. The above expres- 
sion for energy conclusively proves that even 
in the presence of fluctuations, large-N ex- 
pansion gives positive corrections to energy, 
monotonically approaching the exact value 
as one scales up the order. We have checked 
for a range of such higher ordered fluctu- 
ations and have found the previous con- 
clusion sacrosanct. A point of some inter- 
est here would be the variation of such an 
approximated energy with respect to the 
strength uo of the anharmonic oscillation for 
a fixed dimension, N=3 say. Fig. 1 shows 
this variation and evidently tells us that 
there is a minimum in the curve much as 
we would expect it to be. The minimum 
also signifies the fact that the results of 
the large-N approximation would be best 




Non-dimensionalized oscillator strength w 



FIG. 1: Variation of non-dimensionalized en- 
ergy E as in eq. (|17f) with the oscillator strength 
u) for N = 3. The dotted line in the figure 
shows the minimum around which the large-N 
approximation gives the best result. 

valid close to the minimum, that is between 
uj = 0.4 — 0.5 as per Fig. 1. 

As a suggestive example, we might look 
at the next higher modification in the poten- 
tial which gives rise to the following cubic 
equation 



4e 4 Ro 3 +3Vke 3 Ro 2 +(2e2-kuj 2 )R +Vke 1 = 

(18) 

Once again, the above equation can be 
solved analytically using Cardan's method 
and it is rather an easy algebraic exercise 
to show that the energy corrections are still 



positive. 

To conclude, we have shown using a per- 
turbed anharmonic oscillator potential that 
a large-N expansion method provides an 
effective approximation scheme in tackling 
quantum mechanical problems. This is ev- 
ident, since the order of corrections as sug- 
gested by this method monotonically con- 
verges towards the semi-classical limit as 
N —>■ oo. The results offer favorable com- 



parisons with numerical and experimental 
data and might be used in more complicated 
quantum many body problems Q| involv- 
ing exact interaction potentials. 
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